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ABSTRACT 



This thesis examines the use of heat for producing clinical lesions 
in tissue and presents the design and analysis of a resistively heated 
surgical probe. The probe temperature is accurately maintained and 
measured by using a Wheatstone bridge. The probe was embedded in an 
agar-water test medium and the temperature field emanating from the 
probe was measured using liquid crystal tapes. The experimental results 
compare within approximately 10% of a two-dimensional numerical solution. 
A one-dimensional theoretical model is developed and results are 
presented to show the relative effects of metabolism and blood flow on 
the temperature field. 
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NOTATION 



Symbol 


Description 


Typical Units 


a 


Tissue thermal diffusivity 


2 

cm /sec 


a' 


Temperature coefficient of resistance 


ohm/ohm/°C 


c 


Tissue specific heat 


cal/gm/°C 


c b 


Blood specific heat 


cal/gm/°C 


I 


Current 


amps 


k 


Tissue thermal conductivity 


cal /cm/°C/sec 


m b 


Volumetric blood flowrate 


gm/cm 3 /sec 


Q 


Heating rate 


watts 


r 


Radial location in tissue 


cm 


r o 


Probe radius 


cm 


R 


Resi stance 


ohms 


S 

m 


Matabolic heat generation rate 


3 

cal /cm /sec 


t 


Time 


sec 


T 


Tissue temperature 


°C 


T b 


Arterial temperature 


°C 


T c 


Critical temperature 


°C 


T o 


Initial tissue temperature 


°C 


T P 


Probe temperature 


°C 


V 


Vol tage 


volts 


P 


Tissue density 


gm/cm 3 


P’ 


Resistivity 


ohm-cm 



Non-dimensional Groups 




0 






T 




3 = 



m b c b r o 
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I. INTRODUCTION 



A variety of surgical techniques are presently used to produce 
clinical lesions. In addition to the familiar mechanical or "scalpel" 
method, surgical procedures employing chemical, ultrasonic, focused 
x-rays, radiation, and cooling and heating agencies are used to destroy 
descrete regions of tissue, normally nervous tissue. The choice of 
methods depends largely on the nature of the clinical problem. Ideally, 
the surgeon desires a technique that possesses the following character- 
istics: 1) Safety, 2) Reversibility, 3) Reproducibility, 4) Sharp 
delimitation, 5) Hemostasis, 6) Flexibility, 7) Simplicity, and 8) 

Speed of application. At times, it is not possible to meet all of the 
above surgical ideals. In such cases, the choice of procedures is 
based largely on an optimization of the ideals, with safety being of 
paramount importance. Gengler [10] presents a summary of the various 
methods, complete with the advantages and disadvantages of each. 

In this thesis, the use of heat for producing lesions is investigated. 
In particular, the design of a resistively heated surgical electrode, 
henceforth referred to as a probe, is presented. Tissue destruction is 
accomplished by using the probe as a heating element to raise the 
temperature of a predetermined region of tissue above a critical level, 
normally taken to be 55°C [2]. 

The unique feature of this probe is that its surface temperature 
can be accurately measured and controlled by using a Wheatstone bridge. 

As a result, the temperature field in the medium surrounding the probe 
can be analytically predicted by solving the heat equation prior to 
applying the probe. This is quite significant when one considers that 
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the volume of tissue destroyed depends solely on the local tissue 
temperature. 

Temperature fields determined experimentally with the probe 
embedded in a 0.30% agar - 99.7% water test medium compare within an 
accuracy of 13% of those predicted using the analytical model. A novel, 
relatively new, temperature sensing device was used in the experimental 
phase of the study. Liquid Crystals [12], a material that changes color 
over a known and well defined temperature range, were employed to 
obtain a two-dimensional, visual display of the transient temperature 
field which developed around the probe. By using two separate orien- 
tations of the liquid crystals, a three-dimensional picture of the 
temperature field was inferred. In addition to the highly desirable 
visual aspect of the liquid crystal material, it also produces a minimum 
disturbance in the test medium. 

The effects of blood flow and metabolism on the size of the lesions 
produced were also investigated using an analytical model. 
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II. SUMMARY OF HEATING TECHNIQUES 



A. DIRECT CURRENT (dc) METHOD 

Localized lesions have been produced in tissue by applying a direct 
current between two electrodes. An active electrode is inserted into a 
predetermined location in the tissue, and an inactive electrode, or 
ground, is placed on some other part of the body. The resistivity of 
the tissue causes electrical energy to be converted into heat, thus 
increasing the local tissue temperature. When the tissue temperature 
is raised above 55°C [2], the tissue is destroyed. 

Although small and discrete lesions can be produced by the dc 
method, it has several disadvantages. Direct current follows the path 
of least resistivity and tends to have preferred paths in tissue, thus 
causing irregular lesions. Gas bubbles, a result of electrolysis , form 
around the active electrode and block the flow of current. In addition, 
dc devices suffer from polarization which can produce transients that 
stimulate muscle and nervous tissue [23]. 

Nerve stimulation does have one advantage. The probe placement 
can be checked by observing the patient's reaction to the stimulating 
effect of the current. However, other nerves, as well as those in the 
lesion, could also be stimulated. 

Herrero [14] states another disadvantage of the dc method: 

"Animals with dc lesions usually had extensive scar tissue both around 
and infiltrated with the tissue". 

It should also be noted that the dc method depends on both the 
thermal and electrical properties of the tissue. The thermal properties 
are known within an accuracy of approximately 5% [7] and are relatively 
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constant throughout the tissue. However, the electrical properties are 
not well known, are anisotropic and tend to vary with location. 

One advantage of the dc method is that the electrodes are small; 
approximately two millimeters in diameter. Therefore the probes can 
be painlessly inserted and removed, and cause little damage to the 
surrounding tissue. 

In summary, although the dc method can produce localized lesions, 
the irregularity of the lesions, control difficulties, and possible 
side effects of the current, greatly diminish its usefulness. 

B. RADIO-FREQUENCY (rf) CURRENT METHOD 

The radio-frequency technique produces lesions in much the same way 
as the direct current method. A long, thin, needle-like probe emits 
low power continuous wave radio-frequency current. An indifferent 
electrode is placed on some other part of the body and acts as a ground. 
As current is passed between the electrodes, the resistivity of the 
tissue causes electrical energy to be converted to tissue internal 
energy, thus increasing the local tissue temperature. As in the dc 
method, lesions are a function of the tissue's thermal and electrical 
properties . 

Gengler [10] presents a detailed analysis of the rf technique. 

In his summary of the method he states: "The radio-frequency method 

allows controlled and predictable lesions to be produced within safe 
limits for physiological and clinical applications." 

Unlike direct current, rf current has no preferred path in tissue 
and no stimulating effect. However, low-frequency alternating current 
does have a stimulating effect in tissue. Therefore, the probe location 
can be checked by passing low-frequency current between the electrodes 
and observing the patient's response. 
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Brodkey and others [3] discuss another v/ay of assuring proper 
probe location. They found that a temperature rangeof 40-49°C reversibly 
blocks nervous function. Thus, after positioning the probe, and raising 
the local tissue temperature to the reversible range, the patient's 
response can indicate if the probe is positioned correctly. 

Like the dc electrode, the radio-frequency probe is small, can be 
painlessly inserted and removed, and causes little damage to the 
surrounding tissue. 

The rf method does have several disadvantages. The local tissue 
temperature may exceed 100°C causing boiling, gas formation, carboniza- 
tion of the tissue, and adherence of coaguluum to the probe [2]. Also, 
defects in the probe insulation have been known to cause irregular 
lesions along the probe tract. However, these difficulties can be 
avoided by carefully controlling the various physical parameters used 
in the rf technique. 

In summary, the radio-frequency method can be used to produce well- 
circumscri bed lesions and the probe location can be checked by two 
methods. However, as in the dc method, lesions are still a function of 
both the thermal and electrical properties of the tissue. 

C. HIGH RESISTANCE HEATING PROBE METHOD 

In this method, a high resistance electrode is inserted into the 
tissue. Current is passed through the electrode and the resistance of 
the electrode causes the probe temperature to increase. Heat is 
conducted av/ay from the probe into the surrounding tissue. A region 
of destroyed tissue is formed when the local tissue temperature is 
raised above approximately 55°C [2]. 
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Theoretically, there is less uncertainty in producing lesions 
purely by heating than there is in the direct current or radio-frequency 
techniques. The dc method produces both thermal and electrolytic tissue 
destruction. As the lesion forms, it offers an increased resistance to 
the passage of current. Thus, the parameters are changing as the lesion 
forms. The rf method minimizes electrolysis but still depends on the 
electrical and thermal properties of the tissue. On the other hand, the 
high resistance probe depends only on the tissue's thermal properties. 

Carpenter and Whittier [4] attempted to produce lesions in animals 
by using a heating electrode. Although their results were poor, they 
concluded that it was due to the crude instruments employed. They 
stated: "It would seem profitable to explore this technique with a more 

finely perfected thermocautery . " 

Gildenberg [11] conducted studies with a high resistance heating 
electrode and a rf probe and found that the heating probe produced less 
variable lesions than the rf probe. In addition, the rf lesions had 
significantly more lesions along the puncture canal than lesions formed 
by the heating probe. 

As discussed earlier, a temperature range of 40-49°C can reversibly 
block nervous function. This technique can be used to check the 
position of the heating electrode. 

Since lesions formed using the resistance heating method depend 
only on the temperature field surrounding the probe, the surface 
temperature of the probe must be accurately measured and controlled. 

This has been orfe of the difficulties of previous resistance probe 
designs. Watkins [23] used a probe with a thermocouple placed in the 
tip. He found that the actual temperature of the probe could exceed 
the recorded temperature by several degrees if the thermocouple was not 
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positioned carefully. Thermocouples placed in the surrounding tissue 
to record the temperature can distort the temperature field and give 
erroneous readings. Gildenberg [11] inserted a thermocouple one 
millimeter from the heating electrode to record the tissue temperature. 

He determined that controlling the current to maintain a constant 
temperature at the thermocouple was no more accurate than maintaining 
a constant current. 

It should also be noted, that since the lesions produced by the 
heating probe depend only on thermal destruction of the tissue, the 
lesions are smaller than those produced by the direct current or radio- 
frequency current method. 

In summary, a properly designed high resistance heating probe can 
produce small, safe and predictable lesions. The uncertainty is less 
than the radio-frequency or direct current method since tissue destruction 
depends only on the thermal properties of the tissue. 
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III. HIGH RESISTANCE HEATING PROBE 



A. THEORY 

It is well known that the flow of electric current against any 
resistance is accompanied by a dissipation of electrical energy in the 
form of heat. The amount of heat generated is expressed by Joule's 
law: 

Q = IV = I 2 R (1) 

Q = heating rate 
V = vol tage 
I = current 
R = resistance 

The resistance of a conductor at ordinary temperature is a function 
of the length and the cross-section of the conductor, varying directly 
with the length and inversely with the cross-section according to the 
following relationship: 

R = p'f (2) 

R = resistance observed 

l = length of sample 

A = cross-section area of sample 

The proportional ity constant pi is called the resistivity of the 
material. Therefore, all other factors remaining equal, the material 
having the highest resistivity will generate the most heat. 

All metals exhibit some change in resistance with temperature. 

For pure metals, there is a definite and constant relationship between 
the temperature and the resistance of the metal. This characteristic 
is used in electrical -resistance thermometers where temperature is 
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indicated by measuring the change of resistance. Since electrical 
resistance measurements can be made with high accuracy on bridge 
circuits, this type of instrument can give precise temperature 
measurements . 

The theory of the resistance thermometer was applied to the design 
of the high resistance heating probe. From the above discussion, a 
given metal will have a particular resistance at a particular temperature. 
The inverse of this was used in designing the probe. That is, imposing 
a particular resistance on a metal should cause the metal to come to the 
corresponding temperature for that resistance. The resistance can 
easily be controlled by regulating the voltage and current. Thus, the 
temperature of the probe can be measured and controlled accurately 
without the use of any external measuring devices. 

The change in resistance with temperature of a metal is given by 
the following expression: 

R 2 = R 1 [1 + a' (T 2 - T-j )] (3) 

F *1 = resistance at temperature T-j 
R 2 = resistance at temperature T 2 . 

The constant a', is the temperature coefficient of resistivity and 
is expressed as ohms per ohm per degree. It should be noted that, in 
general, a' is a function of temperature. 

B. DESIGN 

Equations (1), (2) and (3) express several factors that were con- 
sidered in designing the high resistance probe. From equations (1) 
and (2) it can be seen that resistivity should be high to generate the 
most heat. Equation (3) shows that the temperature coefficient of 
resistivity should also be high so that the change in resistance with 
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temperature will be large and therefore easier to measure and control. 
Ideally, a' should be constant over the range of temperatures to be 
measured. 

Several metals were investigated. Their properties are listed in 
Table I. In addition, the properties of thermistors were also 
investigated. Although they are small, rugged, can be made in any shape 
and have a long life, their change in resistance with temperature is 
extremely nonlinear. Since this makes temperature control difficult, 
thermistors were not considered further. 



TABLE I 

Resistivity Coefficient 



Metal 


(at 20°C) 


of Resistivity 


Constant 


Platinum 


10 microhm/cm 


.0038 ohm/ohm/°C 




Nickel 


7.8 


.006 




Nichrome 


100 


.0004 




Copper 


1.72 


.0038 




A1 uminum 


2.8 


.0039 




Balco 


19.9 


.0045 


25-1 00°C 


Iron 


10 


.005 





After comparing the properties of the metals in Table I, Balco was 
chosen. Balco is the commercial name of an alloy, consisting of 70% 
nickel and 30% iron, that has the properties that are desirable in a 
heating electrode. Its temperature coefficient of resistivity is large 
(a' = 19.9 microhm/cm) and nearly constant between 25°C and 100°C [24]. 

Electrode insulation was an important consideration in the probe 
design, since a voltage leak in the probe can cause tissue electrolysis. 
Balco is available with a standard enamel insulation. The enamel is 
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made from an oil modified expoxy based resin and has a continuous duty 
operating temperature of 150°C [24]. 

Various shaped heating probes were considered. However, to 
facilitate the construction of the probe, a cylindrical shape was chosen. 
A probe length to diameter ratio (L/D) of ten to one was used to reduce 
end effects and to ensure good experimental data. 

In order to increase the extent of the temperature field in the 
test medium, thus increasing the accuracy of the experimental data, the 
probe was made approximately four times larger than an actual surgical 
probe. 

The heating electrode (see Figures 1 and 2) consists of a skeleton 
of 0.3894 centimeter (cm) diameter glass tubing. Two 0.038 cm copper 
wires were inserted in the tubing. One wire was passed out the tip 
and high resistance 0.0068 cm Balco wire was attached. The Balco wire 
was tightly wound back along the glass tubing for approximately 3.9 cm. 

It was then led back into the tubing through a small hole and attached 
to the second copper wire. The holes at the tip and side of the glass 
tubing were closed with expoxy. The probe was then given a thin 
coating of General Electric No. 7031 Adhesive and Insultating Varnish 
thinned with Xylene to ensure against any insulation defects in the 
Balco wire. The outside diameter of the completed probe was 0.4145 cm. 

C. TEMPERATURE CONTROL 

Precision resistance measurements can be made by using a Wheatstone 
bridge. The bridge (see Figure 3) consists of four resistance arms, a 
voltage source and a detector. When A is to B, as X is to the standard, 
there will be no potential difference bewteen points d and b. Therefore, 
whenever the detector reads zero voltage from point d to b, the bridge 
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Figure 1. Construction of Probe. 



22 







Figure 2. Completed probe 
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Figure 3. Wheatstone Bridge. 
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is in balance and the unknown resistance (X) is obtained by using the 
relationship: 

X = ( £ ) x (Standard) (4) 

The probe was calibrated using a Rosemont variable temperature oil 
bath and commutating bridge. The calibration curve is shown in Figure 
5. Because of the experimental apparatus, a probe temperature between 
45°C and 50°C was desired. The reason for this choice will be discussed 
later. It should be noted, however, that in actual use, a probe 
temperature between 85°C to 95°C would be used. 

As discussed earlier, the flow of current against a resistance 
causes heat to be generated. There is also a change in resistance with 
temperature. Since the temperature of the probe must be controlled 
precisely, it was necessary to ensure that the bridge resistances would 
not change as heat was generated in them. Therefore, 56.2 ohm resistors 
rated at five watts were used for resistances A and B. Referring to 
equation (4), it can be seen that since the ratio of A to B is one, the 
standard resistance will also be the resistance of the unknown (probe). 

The standard resistance was wound using Evanohm wire. Evanohm is 
a commercial wire used in precision resistance standards and is an alloy 
of 75% nickel, 20% chromium, 2.5% aluminum and 2.5% copper. It has a 
low temperature coefficient of resistance over a wide range of 
temperatures [24], The completed Evanohm standard had a resistance of 
356.33 ohms which corresponded to a probe temperature of 48.5°C (see 
Figure 5). 

The final probe circuit (see Figure 4) consisted of the two 56.2 
ohm resistors, the Evanohm standard and the probe. A Lambda Model LP 
413 FM regulated power supply was used to provide power to the probe. 

A Weston Model 1240 digital voltmeter was used as the detector. 
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Figure 4. Temperature Control Circuit. 
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Figure 5. Calibration Curve. 
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Since the probe had a high temperature coefficient of resistance, 
the probe resistance was constantly changing as it generated heat. 
However, by manually regulating the Lambda power supply the bridge could 
be kept balanced and a probe temperature of 48.5°C maintained. 
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IV. EXPERIMENTAL APPARATUS 



A. TEST MEDIUM 

A region of nervous tissue is quite complex. However, to test the 
heating probe, a simplified medium that had the heat transfer character- 
istics of tissue was desired. Since the lesions produced by the heating 
probe are a result of thermal conduction, heat transfer by convection 
had to be negligible in the test medium. In addition, the medium had 
to be clear so that the region surrounding the probe could be seen. 

Grayson [13] suggested that the water content is a major factor 
in determining the thermal conduction of animal tissue. Spells [22] 
investigated thermal conduction in both fluids and tissue and found 
there was a linear relationship between the conductivity and the water 
content of the tissue. The linear relationship was particularly close 
when the water content was above 50%. 

Cooper and Trezek [7] obtained results that indicated the thermal 
properties of human tissue can be predicted within 5% by measuring the 
water content and assuming the remaining portion of the tissue is 
composed of equal amounts of fat and protein. They also found that the 
percentage of water in human tissue varied between 71% and 84%. 

The above results indicate that an ideal test medium should have a 
high percentage of water. Gels would appear to be perfect. Gel 
consisting of gelatin mixed with an appropriate percentage of water 
would approximate tissue as a mixture of protein and water, and thermal 
properties close to those of human tissue could be expected. However, 
a characteristic of gels is that they melt at approximately 22°C. Since 
this would lead to convection effects, gel was not used as the test 
medium. 
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It was found, however, that agar-water solutions consisting of as 
little as 0.25% agar are clear, solidify at room temperature, and remain 
solid up to a temperature of approximately 85°C. In addition, the thermal 
properties of a 1%, or less, agar-water solution are almost identical 
with those of pure water, and within 10% of the thermal properties of 
most human tissues [7]. For these reasons, a mixture of 97.7% water and 
0.3% Difco Bacto-Agar was chosen as the experimental test medium. 

B. LIQUID CRYSTALS AND TEMPERATURE MEASUREMENT 

Under appropriate conditions of temperature and pressure, matter 
can exist in a solid state, isotropic liquid state or gas state. The 
majority of organic and inorganic compounds are comprised of a regular, 
three-dimensional array of molecules or ions in the solid state, i.e., 
a crystal . 

When the crystal lattice is heated, thermal agitation makes it 
impossible for the molecules to maintain a regular arrangement. At a 
given temperature, the crystal begins to melt and changes to a liquid. 

This change, known as the melting point, is usually very sharp and a 
characteristic of the compund. At still higher temperatures, the liquid 
boils and passes into the gas phase. 

Near the turn of the century it was discovered that some compounds 
had a well defined melting point but the melt was opaque. The opacity 
disappeared when the temperature was increased to a definite limit, 
given a true isotropic liquid. Although this intermediate phase was 
fluid to some extent, it also exhibited some characteristics of a 
crystalline solid. This phase was called a mesophase and compounds 
exhibiting this phenomenon became known as liquid crystals [12]. 
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There are three types of mesophase [12]. 

(i) smectic - molecules have a parallel orientation and are in layers 
with their axes at right angles to the layers. 

(ii) namatic - molecules have a parallel orientation of their axes. 

(iii) cholesteric - molecules have a parallel orientation and are in 
layers with their axes in the plane of the layer. 

The cholesteric mesophase is quite different from the smectic and 
nematic and possesses unique optical characteristics. Cholesteric 
liquid crystals having a certain structure exhibit dramatic changes in 
color with very small changes in temperature. This color-temperature 
relationship is constant and reappears everytime the liquid crystal 
passes through a particular temperature range. 

Dowden [9] discusses the structure of cholesteric liquid crystals. 
These crystals are derived from cholesterol or other sterol systems 
and have a basically planar structure. The axial direction in each 
successive layer is displaced slightly from that in the previous layer 

O 

giving a spiral trace. The thickness of each layer is about 3A and 

O 

the pitch of the helix is about 5500A . 

The property which is the basis of color behavior is circular 
dichroism. An incident beam of unpolarized light is split into two 
components having electrical vectors rotated in opposite directions. 

One component is transmitted and the other is scattered. The scattered 
light normally has a wavelength peak that occurs in the visible region 
of the spectrum. 

Temperature changes disturb the intermolecular forces and cause a 
shift in the molecular structure. This changes the wavelength of the 
light that is scattered. Thus, cholesteric liquid crystals exhibit 
dramatic changes in color with very small changes in temperature. 
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Since the colors scattered by liquid crystals are unique for specific 
temperatures, the quantitative measurement of temperature is possible 
to an accuracy of less than 0.1°C [8]. 

Liquid crystals can be applied directly to a surface. However, a 
dark uniform layer is normally placed on the surface before the liquid 
crystals. This provides a standard surface texture, intensifies the 
resulting colors and isolates the crystals from surface contamination. 

Liquid crystals have some application problems. Since they are 
viscous, applied film thicknesses can easily vary and result in non- 
reproducible thermal patterns. Direct exposure of the crystals to 
solvent vapors, oils, dirt and ultra-violet light can cause variations 
in their sensitivity. However, it is possible to encapsulate the liquid 
crystal in a tiny capsule (10 to 30 microns), thus shielding it from 
most of these detrimental effects [17]. 

Encapsulated crystals are available commercially in a variety of 
forms that cover a wide range of temperatures. Rochrome liquid crystal 
tapes manufactured by Hoffman-LaRoche Inc. were used. These crystals 
were coated on a mylar subtratum with an absorptive black background. 

The overall thickness of the crystals and mylar was less than 0.01 
inches. 

The Rochrome crystals are calibrated with an accuracy of 0.1 °C and 
have a response time of less than one second [18]. Parker [19] found 
that the response time of the encapsulated crystals depends on the 
crystal thickness. For a thickness of 0.0035 inches, response time was 
of the order of 200 to 300 milliseconds. 

By placing the liquid crystal tape in the test medium, the tempera- 
ture field around the probe could be seen and temperatures throughout 
the field accurately measured. 
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Rochrome liquid crystal tape lot no. AS-70-36 was used in the actual 
experiments. The color to temperature relationship of this tape was: 

29.9°C - red 

31 .5°C - green 

34.2°C - blue 

During preliminary experiments with the liquid crystal tapes it was 
found that the water in the test medium caused the encapsulated crystals 
to separate from the mylar backing. To prevent this from happening, the 
tapes were given a thin coating of Varathane and the edges were sealed 
with a rubber adhesive. 

C. TEST CELLS 

Two test cells were used to study the three-dimensional temperature 
field surrounding the probe (see Figure 6). A circular container, five 
inches in diameter and four inches deep, was constructed of plexiglass 
and was used to visually observe the radial temperature field. A 
removable ring was fitted into the cell to support the liquid crystal 
tape (see Figure 7). The tape was attached to a thin sheet of mylar 
that was streched across the ring. The probe was inserted through the 
bottom of the cell and was held in place by a bushing that allowed the 
height of the probe to be adjusted. 

The second cell, a 4 5/8 x 5 1/2 x 6 inch rectangular container, 
also constructed of plexiglass, was used to observe the axial temperature 
field near the probe. The probe was inserted through the side of the 
container. A removable bracket was manufactured to hold the liquid 
crystal tape lengthwise along the probe (see Figure 7). Because of 
symmetry, the temperature field on only one side of the probe was 
measured . 
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Figure 6. Test Cells 
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Figure 7. Removable Support Rings for Liouirl Crystal Tape. 



To ensure that the container walls would not disturb the tempera- 
ture field around the probe, both containers were constructed so that 
their dimensions were large compared to the probe dimensions. 
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V. EXPERIMENTAL PROCEDURE 



The test cells were filled with the 0.30% agar-97. 79% water test 
medium for the experimental runs. The probe was positioned in the 
circular test cell so that there was an equal length of probe above and 
below the liquid crystal tape. Thus, the temperature field around the 
midpoint of the probe was measured. The probe was covered by the agar- 
water medium so that the top of the probe was at least five probe 
diameters below the surface. 

All test runs were conducted at room temperature using Rochrome 
liquid crystal tape lot no. AS-70-36 to measure the temperature field 
around the probe. As mentioned earlier, this tape had the following 
temperature to color relationship: 

29.9°C- - red 
31.5°C - green 

34 ,2°C - blue 

A narrow strip of Series 9 Rochrome tape was placed in the test 
medium to measure the medium's ambient temperature. The temperature to 
color relationship of the Series 9 tape was: 

23.0°C - red 

24.0°C - green 

27 .0°C - blue 

The ambient room temperature, normally about 23.5°C, was used as a 
rough check of the temperature obtained from the Series 9 tape. 

Before the containers were filled with the agar-water test medium, 
four reference marks, 4.57 centimeters apart, were placed on the 
corners of the tape.. Figure 8 is a top view of the circular test 
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Figure 8. Top View of Circuit Test Cell. 










cell showing the probe in place, the narrow strip of Series 9 tape, 
and the four reference marks. 

The resistance of the probe was controlled manually by adjusting 
the Lambda power supply to maintain zero voltage on the detector. 
Figure 9 is an overall view of the apparatus. 

Photographs of the temperature field surrounding the probe were 
taken using a Graf lex 4x5 Press camera fitted with a Polaroid Land 
film holder. Polaroid 4x5 Land film Type 57, 3000 speed was used. 

Photographs were taken as various increments of time. The outer 
radius of the temperature field, corresponding to 29.9°C, and the 
reference marks were measured in each picture. Using the actual 
distance between the reference marks (4.57 cm) and the probe radius 
(0.2076 cm), the non-dimensional radius was calculated. 
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Figure 9. Overall View of Apparatus. 






VI. THEORETICAL ANALYSIS 



A. ONE-DIMENSIONAL MODEL 

1 . Metabolic and Blood Flow Effects Neglected 

In order to simulate the heat transfer process in the region 
surrounding the probe, a one-dimensional heat transfer model must be 
developed. In this model it will be assumed that the tissue is an 
infinite homogenous medium with known conductivity, density and specific 
heat and that it is initially at a uniform temperature of T . The 
probe is assumed to be a cylinder of infinite length and radius r Q with 
a constant temperature of Tp (see Figure 10). Thus, the problem is to 
solve for the temperature distribution in a region bounded by a circular 
cyl i nder. 

The general conduction equation is:. 




k = thermal conductivity 

p = density 

c = specific heat 

a = — = thermal diffusivity 

pc 

It is now convenient to introduce the following nondimensional 
variables : 

R = — (distance) 
o 

t = ~~ (time) 
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Figure 10. Mathematical Model of Probe. 
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e 



(temperature) 



T - T o 

= V T o 

Substituting these variables into equation (5) yields, 

30 



1 3 ( R 3r) 



30 _ 

3t 



( 6 ) 



R 3R 
The boundary conditions are: 

(a) at the probe surface (r = a , or R = 1 ) , T = Tp and 8 = 1. 

(b) at large distances from the probe the temperature of the tissue 
remains constant; that is, as r ■+■ T T q or nondimensional ly , 
as R -*■ °° , 0 o . 

The initial condition is at time t = o, T = T q or 8 = 0 for all 
locations. 

Carslaw and Jaeger [5] give the closed form solution of equation 

( 6 ): 



6 = 1 



-- f 

* J 



e xU C Q (u , Ru) du 
u [j l (u) + y2 (u)] 



(7) 



C Q (u, Ru) = J Q (u) Y q (Ru) - Y q (u) J q (Ru) 

where, J Q and Y Q are Bessel functions. Equation (7) was evaluated 
numerically by Jaeger for various values of t and R [15]. 

2. Metabolism and Blood Flow Effects Included 

Pennes [20] was the first to account for metabolic and blood 
flow effects in the heat equation. Pennes suggests that the equation 
take the form: 

3/r 3T\ 
k \ 3 r/ 



3r 



+ m b c b (T b - T) + S m = pc fl 



( 8 ) 
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where : 

= blood flow rate per unit volume 
S m = metabolic heat rate per unit volume 
= blood temperature 

In this model, the effect of blood flow is treated as a perfect 
heat exchanger. That is, it is assumed that the blood enters a unit 
volume of tissue at the arterial temperature, T^, and acquires the 
tissue temperature as it exits the volume. 

If the metabolic heat generation rate is assumed to be constant, it 
can be combined with the blood flow term as: 

m b c b < T o - T) 

where 

S 

T = T. + — 

o b m b c b 

This corresponds to redefining the reference temperature. Cooper 
and Trezek [6] show that for brain tissue, T Q is on the order of 0.5°C 
higher than the arterial blood temperature, T b . 

Equation (8) may be normalized by defining a nondimensional "blood 
flow parameter": 

m c r^ 

6 = — — £ — (blood flow term) 

Substituting this variable along with R, t and e into equation (8) 
yields : 

dt R 90\ 

l l -dr - 66 ■ It < 9 > 

The boundary conditions and initial condition are the same as were 
used for equation (6) . 

Equation (9) was solved with the use of Laplace Transformations and 
the following closed form solution was obtained: 
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2 

TU 



K o (/T) ’ 11 ’ 0 J [u+B/u] [j2 (u) + y2 (u)] 



7T 



(10) 



C Q (u, Ru) = J Q (u) Y Q (uR) - J Q (uR) Y Q (u ) 



Equation (10) was solved numerically for various values of R and t 
and for values of 6 = 0.0, 0.1, 0.5 and 1.0 and the results are tabulated 
in Tables II - XI. The reader is referred to Appendix A for details. 

B. TWO-DIMENSIONAL MODEL 

The one-dimensional model assumed that the probe was infinite in 
length. However, the probe consists of a constant temperature cylinder 
of finite length and a cylindrical insulated stem. The one-dimensional 
solution does not take end effects into account. Thus, it was desired 
to have a two-dimensional solution that would include these effects. 

A digital computer program, by the code name of Trump [16], was 
used to determine the two-dimensional transient temperature field 
around the probe. 

Since TRUMP obtains a solution in dimensional form, a nodal netword 
(see Figure 11) had to be constructed for the geometric configuration 
of the probe. Due to symmetry, it was necessary to generate the field 
for only one half of the cross-sectional area. The dimensions of the 
network were chosen large enough so that the boundary conditions at 
R = 30 and X = 24.0 r Q did not affect the temperature field during the 
transient period which was investigated. The network therefore 
simulated an "infinite" field. 

The TRUMP solution gave values of 0 at each node point for various 
values of nondimensional time. The results are given in Tables XII - 



XVI. 



45 



•30.0R 



a a cr g: g: or a 
Q o ooo o o 0 
oo r< co irW oo c\j 5: 



a 

O 

O' 

c\j 



Cl 

O 



Cl 

O 

O 



, 

222 


\ V r 

00 

OJ 


1 

, V\Z 

T . 


r 

Olif 


20 S 


CMC 

n 

OJ, 


103 

1 


£ 


?! 


•dp j 


00 

r- 

T 


-r 

fA 


3 

n- 


o 


i 

1 


1 

£ ZZ 

1 
















i 










.! 


5 

T 


V 

OJ 

C } 
















- 1 


, 




- * 




, 

| 


po 

C J 


— 1 — 

225 ' 

1 

1 


— T— 

OJ 

OJ 


O 

t“ 

OJ 


00 ■ 

OJ 


o " 

C l 


r m * 
o 

OJ 


201 


rC 

O' 


cr , 1 


l f 

5U 1 

3 

T . 


'.T 

»• 


a* 


>izz_. 


ri 

rv 


cn 

ip 


s 


— r — 

^ ! 


00 

u ^r 

• — 


C\J 

JT 


_ O , 

n 


V_ ‘ 

~ m - 
* — 


■J$i 


.a 


Loj. 

1 1 j 

L 


CL-:. 
: “ o 


0C 


c*lG 


w ......... 

3 
















•j 










rv 


— i — 

62 

L^_ 
















. J 


L ! 








03 

rv 


/ 

* 

* 

/ 

* 


CO 


















| 1 

i- « 

j 


. ■ 








.cn 

r- 


T 

JO 
















■ H 


( 






b 

tP 


Q 

iSl 

L 

£5 


m 

CD 


cn 

lO 

r ■ 


cn 

in H 

T 


NT ' 

T 


‘ cf“ 

v— 


2 1 


c*\ 

of 

4~“. 

c£ 

r" 


I s - 




In 

Cj 


c 


rv 

cn 




CO 

cy 


CO 


■ CO 

in 


in 


" CO ’ 




^ co 


j 1 
£ r L" 

C\|OJ 


CO 


- 


CO 




1 










• ■ 










[_ 








CM 


cn 




















L . 








TO 






- 




- 












1 








V 


to 

n 

1 \ v 


G 

h 

\ \ 


65 


. o 

CD 


in 

in 


. o 

lf> 


. in. 


*1 


r 


j 1 

rt ; 


O' 

AJ 


n 

T - 


- 


n 



Figure 11. TRUMP Nodal Network. 



46 . 



23.0 21.0 19-0 17.0 15.0 12.5 10.5 So 7.0 5.0 3-0 1.0 G 



VII. RESULTS 



The preceding section gave details of the three models that were 
used; a one and two-dimensional model without blood flow and metabolic 
effects and a one-dimensional model including these effects. Numerical 
results for the one-dimensional models are given in Tables I I -XI . 

Tables XII - XVI give the numerical results for the two-dimensional 
model. Values for e are given at various nondimensional locations (R) 
and times (r). The experimental results, differences between the one 
and two-dimensional models, and the relative effects of blood flow and 
time will now be discussed. 



A. EXPERIMENTAL RESULTS 



Figures 12, 13, 14, 15 and 16 are representative pictures of the 
radial temperature field surrounding the probe taken at t = 0.0, 1.0, 
3.0, 5.0 and 10.0. Figure 17, 18, 19, 20 and 21 depict the axial 
temperature field at r = 0.0, 1.0, 3.0, 5.0 and 10.0. The latter 
pictures also show end effects. 

It was convenient to use the following nondimensional variables in 
calculating the experimental results: 



T - T 

e c = y — t~Y~ (critical temperature) 
where, T is the temperature corresponding to a given liquid crystal 

V 

color. In the experimental runs, this color was red and corresponded to 

a temperature of 29.9°C. In a surgical application, T c is the 

temperature at which tissue destruction begins. 

r 

R = — (radius of lesion) 
c r o 

The experimental data and results are given in Table XVII. 
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Figure 12. Radial Temperature Field for 
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Figure 13. Radial Temperature Field for t 
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Figure 14. 



Radial Temperature Field for t 
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Figure 15. 



Radial Temperature Field for x 



5.0. 
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Figure 16. Radial Temperature Field for t 



10 . 0 . 
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Figure 20. 



Axial Temperature Field for t 



5.0. 
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Figure 21 . 



Axial Temperature Field for t 



10 . 0 . 
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The experimental results were compared with the two-dimensional 
model. Figure 22 is a typical graph of the nondimensional time for a 
critical temperature, e c = 0.200. The solid line represents the two- 
dimensional solution and the dashed lines represent the error bound. 

The reader is referred to Appendix B for details of the error analysis. 

The experimental results fall within the error bounds except at 
short times. This is due to the fact that the analytical solutions 
were obtained assuming a step input. That is, it was assumed that the 
probe came to its final temperature immediately. Because of equipment 
limitations this was impossible to duplicate in the experimental runs 
and it took the probe approximately 15 seconds to reach its final 
temperature. 

Figure 23 is a graph of the nondimensional lesion radius versus 
nondimensional time for a critical temperature of e c = 0.235. The 
solid line represents the two-dimensional solution. 

The axial temperature field is shown in Figure 24. Again, the 
solid lines represent the two-dimensional solution. Only one experi- 
mental run was conducted using this configuration. Although the 
experimental isotherms exhibit the same shape as the analytical iso- 
therms, they are lagging in time. A great deal of difficulty was 
encountered in placing the liquid crystal sheet in a uniform plane next 
to the probe and, in fact, the edges of the sheet tended to sag after 
a short period in the agar-water solution. More work is needed on 
this particular aspect of the experiment since experimental errors are 
obviously large. 
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Figure 22. R vs T for 0 = 0.200 and 0 
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Figure 23. R vs t for e„ = 0.235 and p 
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t and e c = 0.285. 
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B. ONE-DIMENSIONAL MODEL 



1 . Effect of Blood Flow on the Temperature Field 

Figure 25 is a graph of the nondimensional temperature versus 
the radial location in the tissue for a nondimensional time equal to 
five (corresponds to a real time of approximately 150 seconds for the 
experimental probe). Four curves are presented: the dashed line 
represents no blood flow effect and the solid lines represent 3's of 
0 . 1 , 0.5 and 1.0. Note from the graph that blood flow tends to decrease 
the temperature field, resulting in a slightly smaller lesion. To gain 
a feeling for the magnitude of this effect, a typical value of 3 will 
be calculated. A typical value for m b c^ of blood flow in the brain is 
0.01 cal/cm /sec/°C [21] and a typical value of the thermal conductivity 
of brain tissue is 0.0014 cal/cm/°C/sec [7]. Therefore, for brain 
tissue, g is approximately equal to 7 r^, where r Q is the probe radius 
expressed in centimeters. The experimental heating probe had a radius of 
approximately 0.2 cm; thus, 3 would be of the order of 0.28. However, 
an actual probe would be approximately one fourth the size of the experi- 
mental probe and 3 would then be much less than 0.1. Figure 26 shows 
that values for blood flow in this range decreased the temperature only 
slightly. 

2. Effect of Increased Time on Temperature Field 

Figure 26 is a graph of the nondimensional temperature versus 
nondimensional time. One can see that at short times, the rate of 
lesion growth is large, but as time increases, the lesion growth 
decreases and approaches a steady size. 

Figure 27 illustrates the effect of blood flow at various times. 
In this graph the temperature at a radius of 2.5 is plotted against time 
for 3 = 0.0, 0.1, 0.5 and 1.0. At short times the blood flow effect 
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Figure 25. 
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Figure 26. 0 vs x for Various Values of R and p 
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Figure 27. e vs x for various values of p and R 



seems to be negligible but as time is increased, the effect of blood 
flow becomes more noticeable. Since the probe radius enters the blood 
flow parameter as a squared term, the blood flow effect would become 
an important parameter for a larger probe at longer times. 

C. TWO-DIMENSIONAL MODEL 

Results for the two-dimensional model are tabulated in Tables XII - 
XVI. The results are given for various positions along the probe 
length, where the X=0 position corresponds to the midpoint of the probe. 

The two-dimensional solution was used to compare the experimental 
results with the theoretical results. Figure 24 shows the isotherms 
along the probe length for various nondimensional times. By extending 
this graph it can be seen that the temperature field around a cylindrical 
probe has an elliptical shape. The dimensions of the ellipsoid can be 
changed by varying the length to diameter ratio of the probe. If this 
ratio is made small enough , the shape of the temperature field will 
approach that of a sphere. 

Figure 28 is a graph of the axial position (X) versus the radial 
location in the tissue for a nondimensional time equal to five and e c = 
0.300. The solid lines represent the two-dimensional solution for an 
L/D of ten, and the dashed line represents the one-dimensional solution. 
Around the midpoint of the probe, the tv/o solutions are in close 
agreement. However, as the distance from the mid-point increases, the 
one-dimensional solution is no longer accurate. 
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Figure 28. Axial Position (X) vs R for x = 5.0 
and e = 0.300. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Previous experimental methods used to predict lesion size and shape 
have generally reduced to trial and error. The usual procedure has been 
to produce lesions in animals, vary the time and current, and analyze 
the results.- 

The theoretical predictions of the temperature field obtained from 
the two-dimensional heat transfer model of the probe fall well within 
the error bounds of the experimental results. Thus, it appears possible 
to analytically predict the size of lesions produced with the constant 
temperature resistance probe. 

Experimental isotherms are reproducible and exhibit sharp delimita- 
tion, particularly at short times. .The probe itself is safe, easy to 

0 

construct and simple. The temperature control circuit is small , 
accurate and easy to operate. 

All results are presented in noridimensional form. However, it 
should be noted that both the experimental and theoretical models are 
simplified and the results have been verified only over a narrow 
temperature range and only for a cylindrical probe with a length to 
diameter ratio of ten. 

The blood flow effect is found to be an important parameter in 
determining the temperature field in tissue if a large probe is used. 

For very small probes, the effect is marginal. Including the effect of 
blood tends to slightly decrease the temperature throughout the field 
and a steady state or equilibrium is reached at an earlier time than if 
blood flow is neglected. As time is increased, the effect of blood 
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flow becomes more predominant. Thus, if heat is applied for any sub- 
stantial period of time, the effect of blood flow is felt more strongly 

B. RECOMMENDATIONS 

It is recommended that additional experimental work be carried out 
with the heating probe over higher temperature ranges and with various 
probe dimensions and L/D ratios. 

As mentioned previously, problems were encountered in attempting 
to place the liquid crystal tape uniformly along the axis of the probe. 
Accurate axial isotherms were therefore not obtained in the present set 
of experiments. Further investigation is needed on this aspect of the 
experiment. 

It is also recommended that a technique for automatically control 1 
ing the probe surface temperature be developed. Two highly desirable 
features of such a control system would be: 1) an ability to select a 
range of probe surface temperatures, and 2) an ability to pre-select 
the duration of the heating period. These two features, together with 
the analytical model, will allow a more accurate prediction of lesion 
size. 
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APPENDIX A: ONE-DIMENSIONAL MODEL SOLUTION 



Statement of Problem: 

A cylindrical probe of radius r Q is embedded in tissue which is 
initially at a uniform temperature T . Within the tissue there is a 
volumetric heat generation due to metabolism of strength $ . In 
addition, there is blood flow in the medium with a flowrate m^, 
specific heat c^ and temperature T^. Assume that the length to dia- 
meter ratio of the probe is large enough so that end effects may be 
neglected and assume the tissue is homogeneous and infinite in extent. 
The probe is at a constant temperature T . Also, assume that at an 

r 

infinite distance from the probe, no temperature change occurs. Solve 
for the tissue temperature assuming constant properties and 
one dimensional radial heat flow. 



Notation : 



Symbol 


Descri pti on 


Typical Units 


a 


Tissue thermal diffusivity 


2 

cm /sec 


c 


Tissue specific heat 


cal /gm/°C 


c b 


Blood specific heat 


cal /gm/°C 


k 


Thermal conductivity of tissue 


cal /cm/sec/°C 


m b 


Volumetric blood flowrate 


gm/cm^/sec 


r 


Radial location in tissue 


cm 


r o 


Probe radius 


cm 


S 

m 


Metabolic heat generation rate 


• ’ 3 

cal /cm /sec 


t 


Time 


sec 


T 


Tissue temperature 


°C 
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Symbol Descri ption 

Arterial temperature 
T 0 Initial tissue temperature 

Tp Probe temperature 



Typical Units 

°C 

°C 

°C 



Nondimensional Groups: 

T - T 

e = j — _ ■ j (temperature) 

P " o 

R = (distance) 

r o 

x = (time) 

r o 

m q ^2 

3 = J L - b . -° (blood flow) 
k 

S 

T 0 = T b + - — — (initial temperature) 
b b 



The general heat conduction equation in cylindrical coordinates is 



3/ r 3T_\ 
k \ 2_r_l 
r 3R 



+ m b c b < T b ‘ T) + S m 



Substituting the above nondimensional 
(A.l) yields: 



pcfl (A.l) 

quantities into equation 



3/ R 30_ 

1 l 3R 

R 3R p u 3t v'-'-'-y 

Equation (A. 2) must be solved in order to describe the temperature 
of the tissue as a function of time. In order to solve the problem, 
boundary conditions and initial conditions must be stated. They are 
as follows: 
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(a) at r = r Q , (at the probe surface), T = Tp or, in nondimensional 
form at R = 1 , 0 = 1 , 

(b) as r -> os, the medium remains unaffected as r -*» ~, T -*■ T Q or, 
in nondimensional form as R -*■ <=», e -*• o, 

(c) at t = o, T = T q or, in nondimensional form at x = o, e = o. 

By using Laplace Transformations, a closed form solution to 

equation (A. 2) can be obtained: 



2 

K 0 m 2 /.» e ' C 0 ( u > Ru > 

6 * — f f al — 5 ? du (A 

K o ^ J ( u + f) [Y o (“) + Y o 

C = J (u) Y (uR) - (uR) Y (u) 
o o o o o 

If the blood flow term is neglected, equation (A. 3) reduces to: 



6 = 1 



2 r ® 

77 oJ u 



" tU C 0 (u, Ru) du 

[ J o («) + (ioi 



(A. 4) 



Equation (A. 4) is the closed form solution given by Carslaw and 
Jaeger [5]. This integral has been solved by Jaeger [15]. 



Solution : 

A finite increment technique was used to solve for the temperature 
field. Equation (A. 3) was rewritten in the following form: 



_ K o (go 
" K 0 (■'T) 




(A. 5) 



e'™ 2 C 0 (u, Ru) 

q (u) = , t : — 5 “ 

(u + £)[J R (u) + y2 (u)] 

Assume an approximation to the integral can be derived by passing 
an interpolating function through several adjacent points and then 
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integrating the interpolating function. Simpson's Rule uses the 

parabola as an interpolatory function between u._-|, u^ and u . + i where 

the u. are evenly spaced values of the interval (u , u ) with a 
i o n 

subinterval length h = u i+1 - u^ . Integrating the parabolic 
interpolating function gives: 




where Y. = q(u. ) 

Expanding equation (A. 6) over the entire interval yields, 

/ q(u) du « | [Y 0 + 4Y, + 2Y 2 + ... 4Y n-1 + Y n ] (A. 7) 

Computer programs 1, 2 and 3 were used to solve this problem for 
$ = 0.0. Programs 4, 5, and 6 were used to solve for 3 = 0.1. To 
solve for other values of 3 only one card must be changed. 

Because q(u) is singular at u = 0, the lower limit of the integral 
was taken as 0.001. The error caused by this was assumed to be negli- 
gible. By varying the upper limit of the integral between 1.0 and 100.0, 
it was found that the solution coverged at approximately 8.0. Therefore, 
an upper limit of 10.0 was used. 

The Y Bessel functions involving the lower limit of the integral, 
0.001, were solved using a separate program and calling the BESY 
subroutine. This was done to avoid underflow and reduce computer time. 

The included programs solve the temperature field for times 
between 0.1 and 1.0 only. To change the time interval only one card 
must be changed. 

Jaeger's results [15] were used as a check on the accuracy of the 
results obtained from Programs 1, 2 and 3 when 3 = 0.0. The greatest 
difference between the solutions was 0.001. 
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APPENDIX B: ERROR ANALYSIS OF EXPERIMENTAL RESULTS 



Three nondimensional quantities were calculated in analyzing the 



experimental results. 

T c - T 0 

e c = j — (critical temperature) 



o 



R = — 
c r_ 



(critical radius) 



x = (time) 
r o 

Ambrosius [1] gives the following equations for the resultant 
errors (R) in computations involving measurements: 

(a) addition 

(A ± a) + (B ± b) 

R = ± (a 2 + b 2 )2 (B. 1 ) 

(b) subtraction 

(A ± a) = (B ± b) 

R = ± (a 2 + b 2 )2 (B.2) 

(c) multiplication 
(A ± a) (B ± b) 

R = ± [(Ab) 2 + (Ba) 2 ]2 (B.3) 

(d) division 

(B ± b) / (A ± a) 

R . ♦ [ < B ^ fl > 2 ; < b > 2 f ( b . 4 ) 

(A) 2 

where A and B are two measurements and an error of ± a is made in the 
measurement of A and an error of ± b is made in the measurement of B 
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The following errors were estimated for the variables in © c : 

T = ± 0.1°C 
c 

T 0 = ± 0.3°C 

T p = ± 0.2°C 

An examination of 6 c along with equation (B.4) indicates that the 
largest error in e c occurs when T q is the greatest, or when the quantity 
T c - T 0 is the smallest. The highest value of T q was in test run 10. 

The resultant errors will now be estimated using data from run 10. 

T q = 25.2°C ± 0.3°C 

T = 29.9°C ± 0.1°C 
c 

T p = 48.5°C ± 0. 2°C 

T - T = 4.7°C 
c o 

T - T = 18.6°C 
P o 

Substituting these values into' equation (B.2), the following errors 
are obtained: 

T - T = ± 0.32°C 
c o 

T - T + ± 0. 36°C 
P o 

From equation (B.4) the resultant error in e c is: 

R = ± 0.025 

or, the percentage error in 6 c = 12.5%. This error would be reduced if 
a higher probe temperature (Tp) or a higher critical temperature (Tj 
were used. 

An analysis of the error in computing R c and t is more difficult 
since the variables in these quantities change with time. It was 
estimated that at short times the error in both R c and x could be of the 
order of 10%. However, this error decreases significantly as time is 
increased . 
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TABLE II 



VALUES OF © FOP NOND I M ENS I ONAL TIME (~C ) ANO NONDI MENS I ONAL 
LOCATION(R) »GN E-DIME MSI ONAL MODEL. 

T=o.l —0.9 

R =1 . 0 — 10.0 



(3 = 0.0 



X 


0 . 1 


0.2 


0.3 


0.4 


R 










1 .0 


1.000 


1 .000 


1.000 


1.000 


1.2 


0.601 


0.691 


0.733 


0.758 


1 .4 


0.316 


0.450 


0.519 


0 . 562 


1.6 


0.143 


0.274 


0.352 


0.405 


1.8 


0.055 


0.156 


0.229 


0.282 


2 .0 


0.018 


0.082 


0.142 


0 . 190 


2.2 


0.005 


0.039 


0.083 


0.124 


2.4 


0.001 


0.017 


0.046 


0.077 


2.6 


0.000 


0.007 


0.024 


0.047 


2.8 




0.002 


0.012 


0 . 027 


3.0 




0.001 


0.006 


0.015 


3.2 




0.000 


0. 002 


0.008 


3.4 






0.001 


0.004 


3.6 






0.000 


0.002 


3.8 








0.001 


4.0 








0.000 


5.0 










6.0 










7.0 










8.0 










9.0 










10.0 











0.5 


0.6 


0.7 


00 

• 

o 


0.9 


1.000 


1.000 


1.000 


1.000 


1.000 


0.776 


0.789 


0.800 


0.808 


0.815 


0.592 


0.615 


0 . 634 


0.648 


0.661 


0.443 


0.472 


0.496 


0.515 


0.532 


0.323 


0.356 


0.382 


0.405 


0.424 


0.230 


0.263 


0.290 


0.313 


0.333 


0.159 


0.189 


0.216 


0.239 


0.259 


0.107 


0.134 


0 . 158 


0.179 


0.199 


0.070 


0.092 


0.113 


0.132 


0 . 150 


0.044 


0.062 


0 . 079 


0.096 


0.111 


0.027 


0.040 


0.054 


0.068 


0.081 


0.016 


0.025 


0.036 


0.047 


0.058 


0.009 


0.016 


0 . 024 


0.032 


0.041 


0.005 


0.009 


0.015 


0.021 


0 . 02 ^ 


0.002 


0.005 


0.009 


0.014 


0.019 


0.001 


0.003 


0.006 


0.009 


0.013 


0.000 


0.000 


0.000 


0.001 


0.001 








0.000 


0.000 
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TABLE III 



VALUES OF © FOR NONO I MENS I ONAL TIME(t) AND NOND I MENS I ONAL 
LOCATION! R) , CNE-DI MENS I ONAL MODEL . 



■ C=l . 0 — 9.0 
R=i .0— 10. 0 
( 3 = 0.0 



~c 

D 


1.0 


ro 

• 

o 


3.0 


4 . 0 


5.0 


6.0 


7.0 


8.0 


9.0 


K 

1 .2 


0.821 


0.854 


0 .8 69 


0.879 


0.885 


0.890 


0.894 


0.898 


0.900 


1.4 


0.672 


0 . 732 


0 . 766 


0.777 


0.789 


0.798 


0.805 


0.811 


0.816 


1.6 


0.546 


0.627 


0 . 665 


0.689 


0.706 


0.718 


0.72 3 


0.736 


0 . 743 


1.8 


0.440 


0 .536 


0.583 


0 . 612 


0.633 


0.648 


0.661 


0.671 


0 . 679 


2.0 


0.351 


0.458 


0.511 


0.544 


0.568 


0.586 


0.601 


0.612 


0.632 


2.2 


0.277 


0.389 


0.447 


0.434 


0.510 


0.531 


0.547 


0.560 


0 . 571 


2.4 


0.216 


0.330 


0.390 


0.430 


0.459 


0.431 


0.493 


0.513 


0.525 


2.6 


0.166 


0.278 


0 . 340 


0.332 


0.412 


0.435 


0.454 


0.470 


0.483 


2.8 


0.126 


0.233 


0.296 


0.338 


0.370 


0.394 


0.414 


0.430 


0.444 


3.0 


0.094 


0 . 1°4 


0.256 


0.299 


0.332 


0.357 


0 . 377 


0.394 


0.409 


3.2 


0.069 


0.161 


0.221 


0.264 


0.297 


0.32 2 


0.344 


0.361 


0 . 3 7 7 


3.4 


0.050 


0.132 


0.190 


0.233 


0.265 


0.292 


0.313 


0.331 


0.347 


3.6 


0.036 


0.108 


0 .163 


0.205 


0.237 


0.263 


0.285 


0.303 


0.319 


3.8 


0.025 


0.033 


0.139 


0.179 


0.211 


0.237 


0 . 259 


0 . 277 


0.292 


4.0 


0.017 


0 . 0" 7 1 


0.118 


0.157 


0.188 


0.213 


0.235 


0.253 


0.269 


5.0 


0.002 


0.021 


0.049 


0 . 076 


0.101 


0.122 


0.142 


0. 159 


0.174 


6.0 


0.000 


0.035 


0.018 


0.034 


0.061 


0.067 


0 . 082 


0.096 


0.110 


7.0 




0.001 


0.005 


0.013 


0.023 


0.034 


0.045 


0.056 


0 . 067 


8.0 




0.000 


0.001 


0 . 005 


0.010 


0.016 


0.024 


0.031 


0.039 


9.0 






0.000 


0.001 


0. 004 


0.007 


0 . 012 


0.016 


0.022 


10.0 








0.000 


0.001 


0.003 


0.005 


0.008 


0.012 
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TABLE IV 



VALUES OF© FOR NOND I MENS I ONAL TIME(x) AMD 
NGNDIMENSIONAL LOCATION(R)t ONE-D I MEN S I ONAL 
MODEL. 









X =10. 


0 -100 


.0 






X 

D 


10.0 


20.0 


R= 1 . 0 
£=0.0 

30.0 


-10.0 

40.0 


50.0 


70.0 


100.0 


K 

1.0 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.2 


0.903 


0.916 


0.922 


0.926 


0.929 


0.933 


0.937 


1.4 


0.820 


0.845 


0.857 


0.864 


0.869 


0.877 


0. 884 


1.6 


0.749 


0.783 


0.800 


0.810 


0.8 18 


0.828 


0. 838 


1 .8 


0.687 


0.729 


0.750 


0.763 


0.772 


0.785 


0.797 


2.0 


0.631 


0.681 


0.705 


0.720 


0.731 


0. 746 


0. 760 


2.2 


0.581 


0.637 


0 . 664 


0.682 


0.694 


0.711 


0.728 


2.4 


0.536 


0. 597 


0.628 


0.647 


0.660 


0.679 


0.698 


2.6 


0.494 


0.561 


0.594 


0.615 


0.629 


0.650 


0.670 


2.8 


0.457 


0.5 29 


0.563 


0.535 


0.601 


0.62 3 


0.644 


3.0 


0.422 


0.497 


0.534 


0.557 


0.574 


0.598 


0. 621 


3.2 


0.390 


0.468 


0.507 


0.532 


0.550 


0.575 


0.593 


3.4 


0.360 


0.441 


0.432 


0.503 


0.526 


0.553 


0. 578 


3.6 


0.333 


0.416 


0 .458 


0.485 


0.505 


0. 53 2 


0.558 


3.8 


0.307 


0.392 


0.436 


0.464 


0.484 


0.512 


0.539 


4.0 


0.283 


0.370 


0.415 


0.444 


0.465 


0.494 


0.522 


5.0 


0.188 


0.277 


0.325 


0.35 8 


0.381 


0.414 


0.446 


6.0 


0.122 


0.207 


0.256 


0. 290 


0.315 


0.350 


0.385 


7.0 


0.077 


0.157 


0.201 


0.235 


0.260 


0.297 


0.334 


8.0 


0.047 


0.112 


0.157 


0. 190 


0.216 


0.253 


0.291 


9.0 


. 0.027 


0.081 


0.122 


0.153 


-0.178 


0.215 


0.254 


10.0 


0.015 


0.057 


0.094 


0. 123 


0.147 


0.183 


0.222 


20.0 


0.000 


0.001 


0.003 


0.009 


0.016 


0.031 


0.053 


30.0 




0.000 


0.000 


0.000 


0.001 


0.003 


0.010 


40.0 

50.0 










0.000 


0.00 0 


0.002 
0. 000 
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TABLE V 



VALUES OF 0 FOR MONO T MENS I ONAL TIMECC) AND NONDIMENS I ONAL 
LOCATION R) T QNE-DI ME NS I ONAL MODEL. 



=0. 1 -0.9 

R= 1 . 0 — 10.0 
^ = 0. 1 



X 

D 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


CO 

• 

o 


0.9 


r\ 

1.0 


1.000 


1 .000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.2 


0.599 


0.688 


0.729 


0.754 


0.770 


0.733 


0.793 


0.800 


0.807 


1 .4 


0.314 


0.447 


0.514 


0. 555 


0.584 


0.606 


0.623 


0.637 


0.649 


1.6 


0.14 2 


0.272 


0.348 


0.398 


0.435 


0.463 


0.485 


0. 503 


0.518 


1.8 


0.055 


0.154 


0.225 


0.277 


0.316 


0.347 


0.372 


0.392 


0.410 


2.0 


0.018 


0.081 


0.139 


0.186 


0.224 


0.255 


0.2 81 


0.302 


0.321 


2.2 


0.005 


0.03° 


0.082 


0.121 


0.155 


0.183 


0.203 


0.229 


0.248 


2 .4 


0.001 


0.017 


0.046 


0. 076 


0.104 


0. 129 


0. 151 


0. 171 


0.189 


2.6 


0.000 


0.007 


0.024 


0.045 


0.06 7 


0. 089 


0. 103 


0.126 


0. 142 


2.8 




0.003 


0.012 


0.026 


0.043 


0.059 


0.076 


0.091 


0.105 


3.0 




0.001 


0.006 


0.014 


0.026 


0.039 


0.052 


0.064 


0.077 


3.2 




0.000 


0.002 


0. 008 


0.015 


0.025 


0.035 


0.045 


0.055 


3.4 






0.001 


0.004 


0. 009 


0.015 


0.023 


0.031 


0.039 


3.6 






0 .000 


0.002 


0.005 


0.009 


0.014 


0.020 


0.027 


3.8 








0. 001 


0.003 


0.005 


0.009 


0.013 


0.018 


4.0 








0.000 


0.001 


0.003 


0.005 


0.009 


0.012 


5.0 










0.000 


0.000 


0. 000 


0.001 


0.001 


6.0 
















0.000 


0.000 



7.0 

8.0 

9.0 

10.0 
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TABLE VI 



VALUES OF© FOP NOND I MENS I ONAL TIMECc) AND NONDIMENS IONAL 
LOCATION(R) , ONE-DIMENSIONAL MODEL. 



T = 1.0 — 9.0 
R = 1 . 0 — 10.0 
£ = 0.1 



D ^ 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


7.0 


8.0 


9.0 


K 

1.0 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.2 


0.812 


0.841 


0.853 


0 . 860 


0.864 


0.867 


0 . 869 


0.9 70 


0 . 872 


1.4 


0.658 


0.710 


0.732 


0.744 


0.752 


0.757 


0.761 


0.764 


0.766 


1.6 


0.531 


0.601 


0.631 


0.648 


0.658 


0.666 


0.671 


0.675 


0.678 


1.8 


0.425 


0.508 


0.545 


0.565 


0.579 


0.588 


0 . 594 


0 . 599 


0.603 


2.0 


0.337 


0.429 


0.471 


0.494 


0.510 


0.520 


0.528 


0.533 


0.538 


2.2 


0.265 


0.362 


0.407 


0.433 


0.450 


0.462 


0.470 


0.476 


0.481 


2.4 


0.205 


0.304 


0.351 


0.379 


0.397 


0.410 


0.419 


0.426 


0.432 


2.6 


0.157 


0 . 254 


0 . 303 


0.332 


0.351 


0.365 


0.375 


0.382 


0.388 


2.8 


0.119 


0.211 


0.260 


0 . 291 


0.311 


0.325 


0.335 


0.343 


0.349 


3.0 


0.089 


0.175 


0.224 


0 . 254 


0.275 


0.239 


0 . 300 


0.308 


0.314 


3.2 


0.065 


0.144 


0.191 


0.222 


0.243 


0.257 


0.269 


0.277 


0.283 


3.4 


0.047 


0.118 


0 .163 


0.193 


0.214 


0.229 


0.240 


0.256 




3.6 


0.034 


0.096 


0.139 


0.168 


0.189 


0.204 


0.215 


C .224 


0 . 231 


3.8 


0.024 


0.077 


0.118 


0.146 


0.168 


0 . 181 


0.193 


0.201 


0.208 


4.0 


0.016 


0.062 


0.099 


0 . 126 


0.146 


0.161 


0 . 172 


0.181 


0.188 


5.0 


0.002 


0.019 


0.040 


0.059 


0.075 


0.038 


0.099 


0.106 


0.112 


6.0 


0.000 


0.005 


0.014 


0.026 


0.036 


0.046 


0.054 


0.061 


0.066 


7.0 




0.001 


0.005 


0 . 010 


0.017 


0.023 


0 . 029 


0.034 


0.038 


8.0 




0.000 


0.001 


0.004 


0.007 


0.011 


0.015 


0.018 


0 . 021 


9.0 






0 . 000 


0.001 


0.003 


0 . GO 5 


0.007 


0. 009 


0.012 


10.0 








0.000 


. 0.001 


0.002 


0.003 


0.005 


0 . 006 
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TABLE VII 



VALUES OF© FOR NOND I MENS I ONAL TIME(-c) AND 
NOND I MENS I ONAL LOCATION(R), ONE-DIMENSIONAL 
MODEL. 

x = 1 o . 0 — 1 00 . 0 

R = 1 .0 — 50.0 
£ = 0.1 



H 


10.0 


R 




1 .0 


1 .000 


1.2 


0 . 872 


1.4 


0.768 


1 .6 


0.680 


1.8 


0.605 


2.0 


0.541 


2.2 


0.485 


2.4 


0.436 


2.6 


0.392 


2.8 


0.354 


3.0 


0 . 31 ° 


3.2 


0.288 


3.4 


0.261 


3.6 


0.236 


3.8 


0.214 


4.0 


0. 193 


5 .0 


0.118 


6.0 


0.071 


7.0 


0.042 


8.0 


0.024 


9.0 


0.014 


10.0 


0 . C07 


20.0 


0.000 


30.0 




40.0 




50.0 





20.0 30.0 

1 .000 1.000 
0.876 0.876 
0.774 0.775 
0.689 0.690 
0.616 0.618 

0.553 0.556 
0.499 0.501 
0.451 0.454 
0.409 0.412 
0.371 0.375 

0.338 0.341 
0.308 0.311 
0.281 0.285 
0.257 0.261 
0.235 0.239 

0.215 0.219 
0.140 0.144 
0.092 0.096 
0.061 0.065 
0.040 0.044 
0.026 0.030 

0.017 0.020 

0.000 0.000 



40.0 50.0 

1.000 1.000 
0.876 0.876 
0.775 0.775 
0.690 0.691 
0.618 0.618 

0.556 0.556 
0.502 0.502 
0.454 0.455 
0.413 0.413 
0.375 0 . 3~ r 5 

0.342 0.042 
0.312 0.312 
0.286 0.286 
0.261 0.262 
0.240 0.240 

0.220 0.220 
0.145 0.145 
0.097 0.097 
0.066 0.066 
0.045 0 . 045 
0.031 0.031 

0.021 0.022 
0.001 0.031 
0.000 0.000 



70.0 10C.0 

1.000 1.000 
0.877 0.877 
0.775 0.775 
0.691 0.691 
0.618 0.618 

0.556 0.556 
0.502 0.502 
0.455 0.455 
0.413 0.413 
0.376 0.376 

0.342 0.342 
0.313 0.313 
0.286 0.286 
0.262 0.262 
0.240 0.240 

0.220 0.220 
0. 1 45 0.145 
C.098 0.098 
0.066 0.066 
0.045 0.045 
0.031 0.031 

0.022 0.022 
0.001 0.001 
0.000 0.000 
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TABLE VIII 



VALUES OF© FOR NOND I MENS I ONAL TIME(-c) AND N OND I MENS I ON AL 
LOCATIONS ) .ONE-DI MENS I ONAL MODEL. 



z 


0.1 


R 




1.0 


1.000 


1.2 


0.592 


1 .4 


0.308 


1.6 


0.139 


1.8 


0.053 


2.0 


0.017 


2 .2 


0.005 


2.4 


0.001 


2.6 


0.000 


2.8 




3.0 




3.2 




3.4 




3.6 




3.8 




4.0 




5.0 




6.0 




7.0 




8.0 




9.0 




10.0 





0.2 0.3 

1.000 1.000 
0.676 0.713 
0.433 0.494 
0.261 0.330 
0.147 0.211 

0.076 0.129 
0.037 0.075 
0.016 0.042 
0.007 0.022 
0.002 0.011 

0.001 0.005 
0.000 0.002 
0.001 
0.000 



“L=0. 1 


— 0.9 


R= 1 . 0 
£> = 0.5 


-*10.0 


0.4 


0.5 


1.000 


1.000 


0.735 


0.749 


0.530 


0 . 5 5 d 


0.374 


0.404 


0.256 


0.289 


0. 170 


0.202 


0.1 10 


0.138 


0.068 


0.091 


0. 041 


0.059 


0.023 


0.037 


0. 013 


0.022 


0. 007 


0.01? 


0.003 


0.007 


0.002 


0.004 


0.001 


0.002 


0.000 


0.001 

0.000 



0.6 0.7 

1.000 1.000 
0.759 0.766 
0.572 0.585 
0.427 0.444 
0.314 0.333 

0.227 0.247 
0.161 0.180 
0.112 0.129 
0.076 0.091 
0.050 0.063 

0.033 0.043 
0.021 0.028 
0.013 0.018 
0.008 0.012 
0.004 0.007 

0.002 0.004 

0.000 0.000 



0.8 0.9 

1.000 1.000 
0.772 0.777 
0.596 0.604 
0.457 0.468 
0.348 0.361 

0.263 0.276 
0.196 0.20° 
0.144 0.157 
0.104 0.116 
0.074 0.085 

0.052 0.061 
0.036 3.043 
0.024 0.030 
0.016 3.021 
0.010 0.014 

0.007 0.009 
0.001 0.001 
0.000 0.000 
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TABLE IX 



VALUES OF © FOR NOND IMENS IONAL TIMECC) AND NOND I MENS I ON AL 
LOCATION(R) , ONE-DIMENS IONAL MODEL. 



X = 1 .0 — 9.0 
R = 1 . 0 — 10.0 
£ = 0 . 5 



*c 

D 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


7.0 


8.0 


9.0 


K 

1.0 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


3 . 000 


1.000 


1.000 


1.2 


0.780 


0.797 


0.801 


0.803 


0.804 


0.804 


0 . 8 04 


0.834 


0.8 04 


1 .4 


0.610 


0.640 


0.648 


0.652 


0.653 


0.653 


S . 654 


0.654 


0 . 6 54 


1.6 


0.477 


0.517 


0.528 


0.532 


0.534 


0.535 


0 . 5 36 


0 . 536 


0 . 536 


1.8 


0.371 


0.419 


0.432 


0.438 


0 .440 


0.441 


0.441 


0.442 


0.442 


2.0 


0.287 


0.340 


0.3 55 


0.361 


0.364 


0.365 


O '. 3 66 


0.366 


0.366 


2.2 


0 . 220 


0.276 


0 . 2 C 3 


0.299 


0.302 


0.302 


5.304 


0.304 


0.304 


2.4 


0.167 


0.223 


0.241 


0 . 248 


0.251 


0.253 


0.2 54 


0.2 54 


0.254 


2.6 


0.126 


0.181 


0.199 


0.207 


0.210 


0.211 


0.212 


0.213 


0.213 


2.8 


0.094 


0.146 


0 .165 


0.172 


0.176 


0.177 


li . 178 


0 . 176 


0.179 


3.0 


0.069 


0.118 


C .136 


0.144 


0.147 


0.149 


D . 149 


0.150 


0.150 


3.2 


0.050 


0.094 


0.112 


0.120 


0.123 


0.125 


3.126 


0.126 


0.126 


3.4 


0.036 


0.075 


0 . 0^2 


0.100 


0. 104 


0.105 


0.106 


0.106 


0.107 


3.6 


0.025 


0.060 


0.073 


0.083 


0.087 


0. 089 


0.090 


0 . 0°0 


0.090 


3.8 


0.018 


0.047 


0.063 


0.070 


0 . 0^3 


0.075 


0.076 


0.076 


0.076 


4.0 


0.017 


0.037 


0.051 


0.058 


0.061 


0.063 


0.064 


0.064 


0.065 


5.0 


0.001 


0 . 01 C 


0.018 


0.023 


0.026 


0.027 


0.028 


0.02 8 


0.029 


6.0 


0.000 


0.002 


0.006 


0 . 009 


0.011 


0.012 


0 . 012 


0.012 


0.013 


7.0 




0.000 


0.002 


0.003 


0.004 


0.005 


0.005 


0.006 


0.006 


8.0 






0.000 


0.001 


0.002 


0.002 


0.002 


0.002 


0.003 


9.0 








0.000 


0.001 


0.001 


0.001 


0.001 


0.001 


10.0 










0.000 


0.000 


0.000 


0.000 


0.001 
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TABLE X 



VALUES OF e FOR NOND IMENS I ONAL TIME(*c) AND 
LOCATION(R) , ONE-DIMENSIONAL MODEL. 

“0=0.1 —0.9 
R*1 .0— 10 .0 
£ = 1.0 



X 

R 

1.0 

1.2 

1.4 
1.6 
1.8 

2.0 

2.2 

2.4 
2.6 
2.8 

3.0 
3.2 

3.4 
3.6 
3.8 

4.0 

5.0 

6.0 

7.0 

8.0 
9.0 

10.0 



0.1 0.2 0.3 

1.000 1.000 1.000 
0.583 0.662 0.695 
0.301 0.417 0.471 
0.135 0.248 0.309 
0.052 0.138 0.195 

0.017 0.071 0.118 
0.005 0.034 0.069 
0.001 0.015 0.038 
0.000 0.006 0.020 
0.002 0.010 

0.001 0.005 
0.000 0.002 
0.001 
0.000 



0.4 0.5 0.6 

1.000 1.000 1.000 1 
0.713 0.724 0.732 0 
0.501 0.521 0.534 0 
0.346 0.370 0.388 0 
0.233 0.259 0.278 0 

0.153 0.178 0.197 0 
0.097 0.120 0.137 o 
0. 059 0.078 0.094 0 
0.035 C.050 0.063 0 
0.020 0.031 0.041 0 

0.011 0.019 0.026 0 
0.006 0.011 0.01S 0 
0.003 0.006 0.010 0 
0.001 0.003 0.006 0 
0.001 0.002 0.003 0 

0.000 0.001 0.002 0 
0.000 0.000 0 



NONDIMENSICNAL 



.7 0.8 0.9 

000 1.000 1.000 
737 0.741 0.744 
544 G. 551 0.556 
400 0.409 0.416 
292 0.303 0.311 

211 0.222 0.231 
151 0.162 0.170 
106 0.116 0.125 
074 0.083 0.090 
050 0.058 0.065 

034 0.040 0.046 
022 0.027 0.032 
014 0.018 0.022 
009 0.012 0.015 
005 0.003 0.010 

003 0.005 0.007 
000 0.000 0.001 
0.000 



TABLE XI 



VALUES OP © FOR NOND I MENS I ONAL TIME HD AND NOND IMENS I ONAL 
LOCATIOM(R) , ONE-DIMENSIONAL MODEL. 



T = 1 .0 9.0 

R= I . 0 10.0 

@ = 1.0 



X. 

n 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


7.0 


8.0 


9.0 


K 

1 .0 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 


1.2 


0.746 


0.755 


0 . 7 56 


0. 756 


0.756 


0. 756 


0.757 


0.757 


0.757 


1 .4 


0.560 


0.575 


0.578 


0.579 


0.579 


0.579 


0.579 


0 . 579 


0. 579 


1.6 


0.422 


0.422 


0 . 445 


0 . 446 


0.446 


0.446 


0 • 446 


0.446 


0.446 


1 .8 


0.317 


0.341 


0.345 


0.346 


0.347 


0.347 


0.347 


0.347 


0.347 


2.0 


0.238 


0.264 


0.269 


0.270 


0.270 


0.270 


0.271 


0.271 


0.271 


2.2 


0.177 


0.205 


0.210 


0. 212 


0.212 


0.212 


0.212 


0.212 


0.212 


2.4 


0.132 


0.160 


0.166 


0.167 


0.167 


0.167 


0. 167 


0.167 


0.167 


2.6 


0. 097 


0.124 


0.130 


0.131 


0.131 


0.131 


0.132 


0.132 


0.132 


2.8 


0.071 


0.097 


0. 102 


0.104 


0. 104 


0.104 


0. 104 


0.104 


0.104 


3.0 


0.051 


0.075 


0.081 


0.082 


0.082 


0 . 082 


0.082 


0.033 


0.083 


3.2 


0.036 


0.058 


0.064 


0.065 


0.065 


0.065 


0.065 


0 • 066 


0.066 


3.4 


0.026 


0.045 


0.0 50 


0.052 


0.052 


0.052 


0.052 


0.052 


0. 052 


3.6 


0.018 


0.035 


0.040 


0.041 


0.041 


0.041 


0.042 


0.042 


0.042 


3.8 


0.012 


0.027 


0.031 


0.033 


0.03 3 


0.033 


0.033 


0.033 


0.033 


4.0 


0.008 


0.021 


0.025 


0.026 


0.026 


0.026 


0.026 


0.027 


0.027 


5.0 


0.001 


0.005 


0.007 


0. 008 


0.009 


0.009 


0.009 


0.009 


0. 009 


6 .0 


0.000 


0.001 


0.0 02 


0.003 


0.003 


0.003 


0.003 


0.003 


0.003 


7.0 




0.000 


0.001 


0. 001 


0.001 


0.001 


0.001 


0.001 


0.001 


8 .0 






0.000 


0.000 


0.000 


0.000 


0. 000 


0. 000 


0.000 



9.0 

10.0 



85 



TABLE XII 



VALUE? OF & FOR NONDIMENS I ONAL TIMECcJ AND NONDIMEN- 
S ION AL L DC AT ION ( R ) » TWO-DIMENSIONAL MODEL ( TRUMP 
SOLUTION) 



t : = o . 5 
R = 1 . 0 
0 = 0. 0 
X =1 P 



50.0 

30.0 



X 

R 

1.5 


0 . 5 


1.0 


3.0 


5.0 


10.0 


20.0 


50.0 


100.0 


0.492 


0.583 


0.696 


0.741 


0.775 


0.791 


0 . 803 


0.830 


2.0 


0.216 


0.330 


0.497 


0.568 


0.622 


0.650 


0.671 


0.723 


2.5 


0.082 


0.176 


0.355 


0.440 


0.513 


0.550 


0 . 576 


0.633 


3.0 


0.027 


0.087 


0.249 


0.339 


0.420 


0.464 


0.499 


0.573 


3.5 


0.008 


0. 039 


0.170 


0.258 


0.353 


0.403 


0.440 


0.512 


4.0 


0.002 


0.016 


0 . 113 


0.194 


0.286 


0.344 


0.390 


0.470 


4.5 


0.000 


0 . OOo 


0.074 


0.142 


0.241 


0.305 


0.350 


0.425 


5.0 




0.002 


0.048 


0 . 105 


0.193 


0.264 


0.316 


0 . 392 


6.0 




0.000 


0.017 


0.050 


0.123 


0.202 


0.258 


0.327 


7.0 






0.005 


0 . 022 


0.077 


0.152 


0.214 


0 . 2 7 8 


8.0 






0. 002 


0.010 


0.048 


0.117 


0 . 179 


0.238 


10.0 






0.000 


0.002 


0.017 


0.065 


0.125 


0.174 


14.0 








0.000 


0.002 


0.015 


0 . 055 


0 . 0 Q 2 


20.0 










0.000 


0.001 


0.013 


0.035 


30.0 












O.COO 


0.001 


0.006 



86 



TABLE XIII 



VALUES OF e FOR MONO I MENS I ON AL TIME ("C ) AND NONDIMEN- 
S ION AL LOCATION(R), TWO-DIMENSIONAL MODEL (TRUMP 
SOLUTION) 

R = 1 . 0 30.0 

-c=0.5 50.0 

£= 0.0 

X = 5R 



T 

n 


0.5 


1.0 


3.0 


5.0 


10.0 


20.0 


50.0 


100.0 


K 

1.5 


0.491 


0.582 


0.695 


0.73? 


0.770 


0.736 


0. 796 


0.819 


2.0 


0.216 


0.330 


0.496 


0. 56 6 


0.617 


0.643 


0.661 


0.699 


2.5 


0.083 


0.176 


0.354 


0.438 


0.506 


0.539 


0.563 


0.611 


3.0 


0.028 


0.087 


0.247 


0.337 


0.413 


0.453 


0.486 


0.541 


3.5 


0.008 


0.040 


0.170 


0.256 


0.343 


0.383 


0.424 


0.485 


4.0 


0.002 


0.016 


0.112 


0. 192 


0.280 


0.332 


0.374 


0.437 


4.5 


0.000 


0.008 


0.073 


0.072 


0.141 


0.232 


0.290 


0.332 


5.0 




0.003 


0.047 


0.103 


0.188 


0.251 


0.297 


0.361 


6.0 




0.000 


0.017 


0 . 049 


0.119 


0.138 


0. 239 


0.300 


7.0 






0.005 


0.021 


0. 07 3 


0.141 


0.196 


0.253 


8.0 






0.002 


0.009 


0.045 


0.107 


0. 163 


0.217 


10.0 






0.000 


0.002 


0.015 


0.057 


0.112 


0.158 


14.0 








0.000 


0.001 


0.013 


0. 047 


0.084 


20.0 










0.000 


0.001 


0.011 


0.032 


30.0 












0 . 000 


0.001 


0.006 
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TABLE XIV 



VALUES OF © FOR MONDIMENS IONAL TIME(X) AND NONDIMEN- 
SIONAL LOCATION(R), TWO-DIMENSIONAL MODEL ( TRUMP 
SOLUTION) 

T=0.5 —50.0 
R=l. 0 — 30. 0 
£ = 0.0 
X=8 . 5R 



X 

R 

1.5 


0.5 


1.0 


3.0 


5.0 


10.0 


20.0 


50. 0 


100.0 


0.485 


0.569 


0.664 


0.700 


0.724 


0.736 


0. 745 


0.772 


2.0 


0.211 


0.315 


0.455 


0.512 


0.552 


0.573 


0. 588 


0.623 


2.5 


0.080 


0.165 


0.313 


0. 380 


0.431 


0.456 


0.477 


0.531 


3.0 


0.027 


0. OPO 


0.213 


0.234 


0.345 


0.376 


0.400 


0.450 


3.5 


0.008 


0.036 


0.142 


0.211 


0.274 


0.310 


0.339 


0.403 


4.0 


0.002 


0.015 


0.093 


0. 155 


0.224 


0.264 


0.295 


0.350 


4.5 


0.000 


0.006 


0.060 


0.112 


0.179 


0.223 


0.258 


0.322 


5.0 




0.002 


0.038 


0. 081 


0.145 


0.195 


0.231 


0.286 


6.0 




0.000 


0.013 


C.037 


0. 088 


0. 143 


0.184 


0.241 


7.0 






0.004 


0.015 


0.053 


0.107 


0.150 


0.201 


8.0 






0.001 


0.006 


0.032 


0.031 


0.124 


0.174 


10. 0 






0.000 


0.001 


0.010 


0.044 


0. 085 


0.128 


14.0 








0.000 


0.002 


0.017 


0.042 


0.073 


20.0 










0.000 


0.001 


0.007 


0.027 


30.0 












0.000 


0.000 


0. 004 
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TABLE XV 



VALUES OF e for NOND I MENS I QNAL TIME(x) ANO NONDIMEN- 
S TONAL LOCAT ION ( R ) » TWO-O IMENS ICNAL MODELdRUMP 
SOLUTION) 

x = 0 . 5 — 50.0 
R = 1 .0 — 30.0 
<?= 0.0 
X=1 1 . 5R 



X 

n 


0.5 


1.0 


3.0 


5.0 


10.0 


20.0 


50.0 


100.0 


K 

1.5 


0.010 


0. 029 


0.094 


0.137 


0.183 


0.22 3 


0.241 


0.293 


2.0 


0.006 


0.022 


0.083 


0 . 124 


0.172 


0.214 


0.231 


0.277 


2.5 


0.003 


0.014 


0.066 


0.106 


0.153 


0.193 


0.214 


0.266 


3.0 


0.001 


0.008 


0 . 050 


0. 086 


0.134 


0.178 


0.197 


0.245 


4.0 




0.002 


0.024 


0.051 


0.095 


0 . 142 


0 . 165 


0.212 


4.5 




0.001 


0.016 


0.037 


0.07 5 


0.125 


0.150 


0. 199 


5.0 




0.000 


0.010 


0.027 


0. 064 


0.111 


0 . 138 


0.184 


6.0 






0.004 


0.012 


0.040 


0.084 


0.118 


0 . 160 


7 .0 






0 . G 01 


0 . 005 


0.024 


0.065 


0 . 098 


0.140 


8.0 






0.000 


0 . 002 


0.014 


0. 049 


0 . 083 


0.123 


10.0 








0 . 000 


0.004 


0.027 


0.058 


0.095 


14.0 










0.000 


0.007 


0 . 027 


0.061 


20.0 












0.001 


0.004 


0.022 


30.0 












0.000 


0 . 000 


0.003 
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TABLE XVI 



VALUES CF © FOR NONDIMENSIONAL TIMEfc) AND NONOIMEN- 
SIONAL LGC AT ION ( R) » TWO-DIMENSIONAL MODEL ( TRUMP 
SOLUTION) 

X=o . 5 — 50.0 
R= 1 . 0 — -30 . 0 
£ = 0.0 
X=1 3 . 5R 



X 

D 


O 

• 

H-* 

• 

o 


3.0 


5.0 


10.0 


20.0 


50.0 


100.0 


K 

1.5 


0.000 0.001 


0.013 


0 . 030 


0.060 


0.102 


0.113 


0.166 


2.0 


0.000 


0.012 


0.027 


0.058 


0.101 


0 . 1 16 


0 . 156 


2.5 




0.010 


0.024 


0.053 


0.095 


0.112 


0.158 


3.0 




0.008 


0 . 020 


0.04 8 


0.090 


0 . 103 


0.147 


3.5 




0.006 


0.016 


0.042 


0.083 


0.102 


0 . 146 


4.0 




0.004 


0. 012 


0.036 


0.077 


0.097 


0.136 


4.5 




0.003 


0.009 


0.030 


0.069 


0 . 092 


0 . 133 


5.0 




0.002 


0.006 


0.025 


0.063 


0.087 


0.124 


6.0 




0.001 


0.003 


0.016 


0.050 


0.076 


0.113 


7.0 




0.000 


0.001 


0.010 


0.039 


0.067 


0.101 


8.0 






0.001 


0.006 


0.031 


0.05 8 


0.091 


10.0 






0.000 


0.002 


0.017 


0.042 


0.073 


14.0 








0.001 


0.004 


0.018 


0.041 


20.0 








0.000 


0.003 


0.004 


0.020 


30.0 










0.000 


0.000 


0.003 
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TABLE XVII. EXPERIMENTAL DATA 



RUN NO. 


X 


*p(cm) 


r (cm) 


T (, o 


T c (°o 


T p (°c) 


©c 


refp 

:cm) p 


ref 

(cm) 


R 


i 


1 • © 


O.IZ 


o.4<S 


22.6 


21.7 


48-S 


0.282 


7. CL 


4-17 


2.X3 


lo- - q | 


2.o 


0-85 


0.53 


22.6 


275 


48-5 


o.28Z 


1 ■ 1,2 


417 


2.57 


C V^O-ci i l ^ 


4 . C> 


1. oc 


o.<G3 


22.C 


27.7 


4«S 


0 .Z 82 


l.cz 


4-77 


3 c>3 




G . o 


1. 16 


0.72 


22.6 


27.7 


785 


0 . 282 


7.62 


4-n 


3 41 




8.o 


1.24 


0.78 


a2.c 


27.7 


7S.5 


0 .X 8 2 


7.62 


4.11 


3.16 




1 0.0 


1.3a 


0.83 


2 2.G 


27-7 


48. S' 


0-2&2 


7.C2 


4-17 


3.H1 


























o.£T 


0.7 ( 


0.45 


23.7 


27.7 


48. S' 


0.244 


7.32 


4.51 


2.14 


u - s-nt 


1. S' 


0 . 8 I 


O-Sl 


23.7 


27-7 


48.5 


0-244 


1.3 2 


4-57 


a. 45 


V~ <x ci { <x V 


3.0 


o.7 7 


0 . Go 


23.7 


27 7 


4G5 


0-24 4 


7.32 


4. si 


2.41 




4.o 


l.ol 


o.67 


2.35 


27 .7 


48- S 


O 244 


7.32 


4-57 


3-2o 




S.o 


1-/3 


o.ll 


23.7 


27.7 


48. S' 


0-244 


1.32 


4.57 


3.4o 




6- o 


/•at 


0.15 


2£.7_ 


27.7 


48-5 


0-244 


1.32 


4.57 


3.63 




1.0 


(-'2G 


o.7& 


23.7 


27 7 


480' 


0-244 


1.32 


4-S1 


3.18 
























3 


Z.o 


0.17 


0-62- 


2 . 3.0 


27.7 


48 5 


o.27l 


1.34 


4-57 


2.7 1 


It- 8 - K 


4.o 


1 ./ 1 


0.73 


23.0 


27.7 


48.5 


0-27 l 


1.34 


4-57 


3.50 


Cv'cvd . o-l ) 


6.0 


t.zz 


0.80 


23.0 


27.7 


4&*S 


O 27 \ 


7.34 


457 


3 . 85 


X- IQ. 


8.o 


!.4o 


o.ai 


23.0 


27 7 


48. S 


0-27 1 


1.34 


4S7 


4.2o 




( 0.0 


I. SO 


0.7 3 


23.0 


27. 7 


48.5 


0-271 


7.34 


4-5 7 


4.47 
























4- 


2.0 


O.?? 


0.62- 


23.0 


27.7 


485 


0.27 \ 


7.32 


4-S7 


J4 7 


u-7 -*7i 


4.0 


1./7 


0.13 


23.0 


27.7 


485 


0.27 \ 


7.32 


4-5 7 


-3.5o 


CA- ci * O- l ^ 


&.o 


1-32. 


0-83 


23.o 


27.7 


485 


o-27 l 


7.32 


457 


3.7 | 


X « iO 


S.o 


1.41 


0.88 


23.o 


27.7 


48.5 


0 . 27 V 


7.32 


4- S7 


4-2.3 




1 0.0 


l-So 


0.73 


23.o 


27.7 


48-5 


0-27 1 


1.32 


4-57 


4.41 
























S' 


2.o 


0. 71 


O.S7 


24. 0 


27.1 


4 8. 5 


0-24o 


1.1 1 


4-57 


2.82. 


n '/7 - 1 1 


4.o 


1-12. 


0.12. 


A4.0 


277 


48.5 


0 - 24o 


111 


4.51 


3.4C 


(jr cxcJU c*_0 






















X- A R. 
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TABLE XVII. EXPERIMENTAL DATA (Cont'd) 



RUN NO 


X 


rp(cnri) 


I"' (cm) 


T’c 


T/c) 


T/o 


©c 


(2. 

< 4 -^ 

<L)E 


ref 

(crrO 


R 


G 


2.o 


o.ll 


o.Sl 


23. C 


29.4 


48.5 


0253 


9-U 


4.59 


2. S3 


i v - Zo - q i 


£ . o 


1.2ft 


o 23 


23. G 


214 


48.5 


O 2S3 


9.1 1 


4. 59 


3 .18 




S.o 


13ft 


0.81 


23.C 


21.1 


48.5 


0-253 


9.1 1 


4.S9 


4.2ft 


A- 1R 


t O. O 


14S 


0.15 


23. G 


21- 1 


49-5 


0-253 


9.11 


4.S4 


























7 


2.o 


O.SG 


0.53 


24.2 


21-1 


48.5 


0.235 


C.S3 


4.5 9 


2.91 


\ 

rv) 


4 .o 


l.ol 


0.13 


24.2 


211 


48.5 


0.235 


6.83 


4.51 


3.51 


(jV“cxcl V (X.V ") 


S.o 


1.31 


o.U 


24.2 


21-1 


49-5 


O. 235 


G-S3 


4-51 


4-4 2 


X--1R 


10.0 


1.42 


0.15 


24.2 


_21.1 


485 


o. 235 


£.<83 


4-51 


4-51 
























8 


\ . o 


O.G 1 


o 4 1 


2-42 


211 


48 5 


0.235 


£.83 


4.59 


1.19 


, iv- 2.i -a i 


3.o 


o.Tl 


o £ 


24- 2 


21.1 


49- 5 


0-235 


£.<83 


4.51 


3 ,2o 


(v*&.citc^O 


zs:o 


in 


o.iS 


24.2 


21.1 


48.5 


0-235 


£.<33 


4.51 


3.99 


>s= i R 


1.0 


1.32 


o.ee 


242 


21.1 


40. 5 


0-235 


G 83 


4 51 


4 2G 




1.0 


144 


o. l£ 


24 2 


21.1 


48-5 


0-235 


£ <3 3 


4 51 


4-G2 
























7 


1.5 


0.14 


0.41 


24.2 


21.1 


49.5 


O 23S 


£.83 


4.59 


239 


H - 2l - K 


3.5 


\.o3 


oGl 


24 2 


21.1 


48-5 


o.23 5 


G83 


4 >1 


3. 3 l 


CV'o.dUa-l^ 


s.s 


LIS 


o. 91 


24 2 


21 1 


485 


0.23 5 


£.83 


4 59 


3.8 2 


^ - 1R 


9.5 


1-3 V 


a.88 


24 2 


2.11 


40-5 


O. 235 


£.33 


l.sn 


4-22 
























1 o 


i . o 


on 


0.4 3 


252 


211 


48.5 


0.2o2 


£.63 


4.59 


2 .ol 


/I-22 -HI 


1.5 


o.lG 


o.5l 


2 s Z 


211 


48.5 


o-i-ol 


£-6 3 


4.51 


2.4G 


CwO 


2 .o 


o.84 


o. SG 


25-2 


21 1 


48.5 


O.Zo2 


£ - S 3 


4 51 


2.10 


X-= VR 


3. o 


\ o 2. 


0.66 


252 


21 1 


48-5 


O- 2 .oZ 


£.83 


4.51 


3.2S 




3.0 


A- A \ 


o. a<s. 


2 52 


21.1 


48.5 


e> ZoZ. 


£><33 


4.51 


358 




A-o 


1. c>l 


o.l3 


2 S 2 


21 1 


49.5 


O 2o2 


GQ3 


4 51 


3.5> 




5.0 


l.ZG 


o 84 


25 2 


211 


48.5 


0-^.0 2 . 


6 83 


4 51 


4.o5 




5.5 


1.32 


o.©8 


252 


2 H 


48.5 


o • £*ocL 


' £.<33 


4.51 


4-2£ 




G .o 


1X4 


o.g 3 


252 


21 1 


48-5 


o 2 o 2 


£•<8 3 


4.59 


4.0 l 




1.0 


I.3S 


ol3 


25- 2 


21 . 7 


48.5 


o. Zo 2 


£83 


4.59 


"4 • 4£ 
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TABLE XVII. EXPERIMENTAL DATA (Cont'd) 



RUN NO. 


X 


£(cm) 


P (cm) 


T 0 (< o 




T p (°c) 


©c 


ref D 

(c m) 


ref 

(cm) 


R 


10 (corxO 


7. S' 


M7 


o 77 


25.2 


27.7 


46.5" 


o- 20 2 


g.83 


4.57 


445 




8.o 


1.40 


o.?4 


2S2 


27.7 


7e.tr 


a 202 


€.E3 


4.S7 


4. SO 




9.o 


I.SZ 


f.oa 


as-. 2 


27.7 


485 


o 203 


S.’SS 


4.57 


4.7 ( 




9.o 


l Go 


1.0 1 


2S-.2 


27.7 


4g.iT 


0.20 2 


6.33 


4.S7 


5. 16 




/c . o 


I.S2 


\.oZ 


25.2 


27.7 


4 £5 


o.2o2 


6.83 


457 


4 7 1 
























il 


2.o 


6.71 


0.46 


22. S 


27- 7 


4-8.tr 


02§5 


3.4o 


2.2 \ 


2.23 


(i - 28- 7 1 


4 .& 


0-7/ 


o 54 


2. Z.s 


27- 7 


48.S 


o. 2St5 


3.4o 


2.21 


Z 


(CCMG.0 


G.o 


i.o2. 


o . G G 


22.tr 


24.7 


49. S 


o '28 s 


3.4o 


2 .2 l 


3.16 




So 


l.o7 


0.67 


22-tT 


27-7 


4S.S 


o. 28S- 


3.4^ 


2.21 


3.34- 




/o o 


i.*2 


0.73 


22. S' 


27-7 


48. S 


o-28S 


3.4o 


2.21 


3 . So 
























x - s.sa 


2.o> 


O.&l 


0.40 


22 >5 


27-7 


7S.S- 


o. 29. ts- 


3 4o 


2. XI 


!7 / 




4.0 


o ■ 66 


0.43 


22-tT 


27.7 


48- S 


o.2SS 


3 .4° 


2.21 


2.o7 




& . o 


o.9l 


o.53 


22 or 


27-7 


48-S 


o. 19^ 


3. 40 


2.2» 


2.S4 




S.o 


o.7 | 


O.S 7 


22 .tr 


27.7 


49. S 


0-28 S’ 


3.4o 


2.21 


2.86 




l o o 


0.77 


o.63 


22. t> 


27.7 




0.285 


3. 4o 


2.21 


3,o2 
























X - loR 


2. o 


o .S' I 


O. 3^ 


22 5 


27.7 


48. r 


0-28 5 


3.4o 


2.21 


1.S7 




4 . o 


O.^G 


O.3<o 


22.tr 


27.7 


40.tr 


o- 20S 


3.4o 


2.2V 


\-7tr 




G * o 


O.GG 


o. 43 


z 2.^ 


27-9 


48- S’ 


o 2J8S 


3. A<=> 


2.2V 


2. .07 




S. o 


0-7 | 


o-46 


22 . r 


24.7 


48. S 


o-ZQ5 


3 .4 0 


2.21 


2.23 




lo.o 


o.qG 


0.47 


22 .tr 


27.7 


40- ^ 


0235 


3.4 o> 


22.1 


2.37 
























*•= W.SG. 


2.o 


0.33 


o.Zo 


22 .r 


27.7 




0 2.85 


3.4 0 


224 


l.c>o 




4.0 


0.3 3 


o- 2o 


22-5- 


27.7 


4®^ 


0 285 


3.4o 


22. 1 


lOo 




G.o 


0.33 


o. 2.0 


22. r 


27.7 


48. r 


o 285 


3. 4o 


22.1 


J.Oo 




S.o 


0-3 3 


O. Zo 


22 .tr 


27.7 


40. S 


o 285 


3. 4 0 


22.1 


l.Op 




\0 o 


O. 41 


0.2.6 


22. s 


27-7 


45*3 


0 2&S 


3 .4 0 


22.1 


1 27 
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COMPUTER PROGRAM 



o 

»— i 

o 



o 



CC 

< 

o 



o 

CC 

Li- 



00 

UJ 



-J 


.u 






< 


< 






oc: 








o 


o 






1 X 1 


LU 


00 




h" 


K 


-J 


si 


~z: 


2 


< 


a: 




i 


> 


UJ 


CO 




or 


l — 


wLL 


LL 


UJ 




aco 


O 


1 — 


3 : 


a: 




z 


O 


UOh 


K 




-J 




►— t 




LL 


zxsi 


S~ 


UL 




LL >-h 




O 


O 


Z-J 


-J 




a 


-JO 




QC 


co 


<T H-i (V* 


a" 


LU 


_j • 


ZO 0 UJ 


Hiuo co 


cco 


tfzsca 


• si 


O II II 



CC 

O 

LL. 



00 
UJ 

O h- 

<£ r-n co vf* co »— h co co co c*** 

>v^h*-cococMcMCMvO(MCMO 
r-HOm^f^cocMcot^<)co 
-jr-r^uoOLn^Hh-co^oo 
LUN$-stCOcO<M(M-~ ihO 

• M-<r mm ( 0(0 

| | M I I 

lu M | l I I | l | || || || || H || 

03 || || || || || || || || || — . — — — 

o ?-H cm co m o r- oo o •-* 

r— I (\I CO LPl \j 0 CO O'' *-H pH *—H **H r 1 r-H r— { r~H f“H r-H CM CM 



como 
•OOOh 

ICD^HO 

)OO^vt 

) O' o n o O' 



CMCMxOcOCMO 
CO •— i CO CO CT — h 
OhOCMAO 
coco O' men cm 

00 xO (V s * »-^ CM^ 
co coco cor- r- 



in cm .-Hr- cm r- 
n3 cm co cm r- co 
•-H o co ro oo o 
p~H co •— i cm coin 
mm —tor-m 

Nf^^'O'O'O 



co in co » — < 
cr^r- vo 
^rnocc 

OOf-H^OO 

m cmo co 
O vO'Oin 



^ coco com co co 

11)111 
It II M II II II 



lu lu o o Q. c o o < < cdcc cc cc cc cc cc r£ a: a: cc 
hSJ OHZMhhCLOOOOOOOCOO 
II II || II II II LUUJIC >->->>>->> >->>- 
LUQ<<coa;^ 2 :cca 3 *-H<<<<<<<<<< 



cccccc.cc 

LDLfYo; 

OOOO 

>->>-> 

<<<< 



CC CC CC cc QC cc 
XCCqCCCcCCC 

oooooo 

>>->>->> 

<<<<<< 



comcofocoro 

Mill) 

iUl 11 11 11 11 

(MCOvfinOh' 
CM CM CM CM CM OsJ 

QCCCCCCCCCCC 

(XCLCCCCCCCC 

OOOOOO 

>■>>>-» 

<<<<<< 



COCO CO CO 

II I I 
II II II II 

coo o^c 

CM CM COCO • 

wrw- wwO 

OCCC Cf.eC II 

cc cc cc cc s: 

ooooo 

<<<<H 



cc 

< 

> 

LU 

zc 



CL 

oo 

OH 



O II 

o *-* 

o 

• aio 

O LU O 
II h-O 
SI 3 »— • 
OO 
oo O 
aih-o 
a: 

< 



00 



o o o o o 



94 



START INNER OD LOOP. RADIUS VARIES FROM 1.0 TO 4.0 FOR EACH 
VALUE OF TIME. 

DO 1100 K= 1 » 31 
R=RSUM + 1 .0 
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